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Symbols 


A 

nucleon number 

d 

overlap distance, fm 

EM 

electromagnetic 

Aphotonmax 

maximum photon energy, MeV 

GDR 

giant dipole resonance 

9n 

neutron branching ratio 

Qp 

proton branching ratio 

Apts 

number of integration intervals 

#o.l 

10-percent-charge density radius, fm 

7th 

threshold energy, MeV 

Z 

proton number 

r 

GDR width, MeV 

^EM 

electromagnetic dissociation cross section, mb 
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Introduction 

As the United States space program heads 
toward the era of permanent occupation of the near- 
Earth environment via the habitation of a permanent 
Space Station, the problem of protection of astro- 
nauts from the harmful effects of cosmic radiation 
assumes greater and greater importance. Of spe- 
cial concern are the effects of galactic cosmic rays in 
the form of relativistic nuclei (ref. 1). When a rela- 
tivistic nucleus impinges upon a spacecraft wall, it 
undergoes several reactions, an important one of 
which is fragmentation, whereby the nucleus breaks 
up into smaller pieces which subsequently decay. It 
is these secondary particles which contribute to the 
radiation environment inside a spacecraft. The prob- 
lem of fragmentation of relativistic nuclei is the sub- 
ject of a continuing research effort by the present 
authors and others (ref. 2). An alternative mech- 
anism, which can produce the same sort of radia- 
tion environment as the fragmentation mechanism, 
is that of electromagnetic dissociation by which the 
electromagnetic field of one nucleus excites states in 
another nucleus which subsequently decay. Norbury 
and Townsend have led a study (ref. 3) of the signif- 
icance of this mechanism to galactic heavy ion break 
up and have concluded that it is of major importance 
(i.e., large cross section) for the removal of a few nu- 
cleons when a relativistic nucleus impinges upon a 
spacecraft wall. Consequently, the mechanism of the 
electromagnetic dissociation must be included when 
one is trying to predict the radiation environment in 
the interior of a spacecraft when the exterior envi- 
ronment is composed of galactic heavy ions. 

In the present work, it is assumed that the 
photonuclear cross section is totally dominated by 
the giant dipole resonance (GDR) (refs. 4 through 
9) and, furthermore, that the GDR decays only via 
single nucleon emission. 

Parameterizations 

The theoretical formulations of electromagnetic 
(EM) dissociation cross sections (ref. 3) are certainly 
adequate as they stand. However, because of the 
large number of parameters that need to be input 
“by hand” into the EM code (ref. 3) in an interactive 
session, one cannot use the EM code as an integral 
part of a heavy ion transport code. Thus, parameter- 
izations are provided for the giant dipole resonance 
width T, the 10-percent-charge density radius Rq i , 
the photoneutron and photoproton threshold ener- 
gies 7th and branching ratios g p and g n , the nucleus- 
nucleus overlap distance d, and finally the number 
of integration intervals used in the numerical 


integration of the photon spectrum and photonuclear 
cross section. 

An adequate parameterization of all of these 
quantities means that the only inputs required for 
the EM code are the projectile kinetic energy and 
the mass (A) and charge ( Z ) numbers of the projec- 
tile and target. 

Giant Dipole Resonance Width 

Inspection of table 1 of reference 3 indicates 
that light nuclei have large widths (of the order of 

10 MeV) and heavy nuclei have relatively small 

widths (around 4-5 MeV). However, attempts to pa- 
rameterize T have not been very successful (ref. 5) as 
it can be sensitive to nuclear structure effects. Fortu- 
nately, the EM cross sections are not very sensitive to 
T (ref. 3), and a very simple parameterization seems 
to be adequate. Thus, we take 

10 MeV (A < 50) (1) 

r - 4.5 Mev (A > 50) (2) 

These values are listed for a representative sample of 
some nuclei in table 1, together with the values from 
reference 3. 

Branching Ratios 

Consistent with our assumption that the GDR 
decays only via single nucleon emission*, the neutron 
branching ratio is 

9n = 1 — Qp ( 3 ) 

where it remains to parameterize g p , the proton 
branching ratio. We shall use the one provided by 
Westfall et al. (ref. 6) as 

g p = Min [Z/A, 1.95 exp(-0.0075Z)] (4) 

which denotes the minimum value of the two quan- 
tities in brackets. Again these values are listed in 
table 1. 

Photonuclear Reaction Thresholds 

The cross section is not strongly dependent upon 
the particle reaction thresholds; therefore, we use 
0.001 MeV to avoid obtaining an infinite value in 
the EM codes (note that a value of 0 MeV would 
cause numerical difficulty). However, the true value 
is on the order of 10 MeV (see table 2), and one would 
therefore think that this would be a poor approxima- 
tion, especially as the dominant contribution to the 
virtual photon spectrum occurs at low energy. This 
effect is offset due to the giant dipole cross section 
becoming negligible near the threshold and therefore 



contributing nothing to the total EM cross section. 
As a result of 

T th = 0.001 MeV (5) 

(see table 2), the code now no longer requires nu- 
clear mass excesses as input; this is a considerable 
improvement in simplicity. 

10-Percent-Charge Density Radius 

Nuclear radii are well parameterized by 
1.18 A 1 / 3 fm. Motivated by this result, we investi- 
gated the possibility of adding a constant to this form 
in order to fit the 10-percent-charge density radii. 
(See table 3 of ref. 3.) The best fit was obtained 
with 

Ro.i = 1.1&4 1 / 3 + 0.75 fm (6) 

and the results are listed in table 3, together with the 
results of reference 3 where a variety of complicated 
models were used to calculate Ro.i- As can be seen, 
equation (6) reproduces the 10-percent-charge radii 
very accurately. 

Overlap Distance 

The nucleus-nucleus overlap distance was chosen 
as 

d = 0 fm (7) 

Integration Variables 

The simple trapezoidal rule was used in the nu- 
merical integration. A very reliable result can be ob- 
tained with the number of integration intervals being 

JVpts - 50 (8) 

Furthermore, just as the lower limit in energy is taken 
as 0.001 MeV, inspection of figures 10 through 30 of 
reference 3 indicates that an upper energy limit 

^photonmax = 70 MeV (9) 

is sufficient to cover the entire GDR region. The 

reader is reminded that the relatively small value 
used in equation (9) results from the dynamics of 
the GDR and is not strongly dependent upon the 
incident kinetic energy of the projectile. 

Results 

Equations (1) through (9) summarize all the pa- 
rameterization of the nucleus-nucleus EM dissocia- 
tion cross sections. The only required inputs to the 
new EM code (see listing in the appendix) are the 


projectile kinetic energy and the charge and mass 
numbers of the projectile and target. 

The theoretical (ref. 3) and experimental GDR 
photonuclear reaction cross sections (refs. 5, 10, and 
11), which appear in figures 1 through 13, are pre- 
sented with these parameterizations. The compari- 
son between the present work and that of reference 3 
is quite good except for 58 Ni and 54 Fe. The case of 
54 Fe is particularly bad because the present param- 
eterization of g p leads to poor results for nuclei far 
from stability. We shall return to this point later. 

The total nucleus-nucleus EM dissociation cross 
sections are presented in tables 4 through 8. (Note 
that table 7 presents results for target fragmentation, 
whereas all the other results are for projectile frag- 
mentation. To calculate target fragmentation with 
the code in the appendix, one should simply swap 
the projectile and target.) Also included in these ta- 
bles are the calculations of reference 3. Experimental 
data (refs. 7, 8, and 9) are also included in tables 5, 
6, and 7. 

The agreement between the parameterization of 
<tem an d the data is reasonably good. However, gen- 
erally the parameterization gives a larger result than 
the projectile fragmentation data of tables 5 and 6 
and a smaller result than the target fragmentation 
data of table 7. The reason for this is that the pa- 
rameters were adjusted to give the best overall agree- 
ment with the data. Note however, that the errors 
associated with the data are very large. The errors 
in the data of Heckman and Lindstrom (ref. 7) in ta- 
ble 5 are typically 50 percent and those of Mercier 
et al. (ref. 9) in table 7 are around 25 percent. Thus, 
even though our parameterization agrees with the 
data usually within 20 percent, this very large uncer- 
tainty in the data means that the parameterizations 
are really as uncertain as the data. Therefore, there 
is a very urgent need for more accurate data. 

The 18 0 (ref. 8) in table 6 has very small errors, 
and our parameterization typically overestimates the 
data. There are two reasons for this. First, the 
photonuclear cross section (figs. 10 and 11) displays 
features very different from the normal giant dipole 
resonance. Thus, any attempt to fit the standard 
GDR is bound to fail. Second, the parameterization 
of the branching ratio in equation (4) is most accu- 
rate for stable nuclei. For unstable nuclei like 18 O 
and 54 Fe (refer to earlier discussion), its accuracy is 
limited. Fortunately, however, the production cross 
section of these rare nuclei is reasonably small, and 
the corresponding inaccuracy will not manifest itself 
as a large inaccuracy in a transport code. Neverthe- 
less, a good description of the electromagnetic cross 
sections for unstable nuclei would be valuable. 
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Finally, as seen from tables 4 through 8, agree- 
ment between our parameterizations and the results 
of reference 3 is excellent. 

Concluding Remarks and Future Needs 

The parameterizations presented herein are able 
to match the experimental data to within the uncer- 
tainty of the data. This work has demonstrated that 
the following points need to be further addressed (in 
order of importance): 

1. Much more accurate data are needed 

2. A more accurate theoretical model is needed 
to describe the photonuclear reaction cross 
sections; this would consequently improve the 
parameterizations 

3. A better description of unstable nuclei (such 
as 18 0 and 54 Fe) is needed, particularly 


the branching ratios and photonuclear cross 
sections 

4. Multiple nucleon emission should be included 

5. Alternatives to the Weizsacker- Williams meth- 
od for obtaining the photon spectrum should 
be investigated 

6. The size of the interference between the strong 
and electromagnetic forces needs to be deter- 
mined (we have assumed that it is zero (refs. 8 
and 12)) 

7. Multipolarities other than electric dipole need 
to be included (refs. 12 and 13) 

8. Curvilinear, rather than straight-line, trajec- 
tories should be considered (refs. 12 and 13) 

NASA Langley Research Center 
Hampton, Virginia 23665-5225 
May 4, 1988 
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Appendix 


Computer Code 

The computer program of reference 3 has been modified to include the parameterizations of the present 
work. The only required inputs are the projectile kinetic energy and the projectile and target mass and charge 
numbers. At the end of the program, a sample output is provided. 

Program Listing 

10 REM COULOMB PARAMETERIZATION 

20 REM 

30 REM 

40 REM 

50 FIXED 2 

60 REM 

70 REM 

80 DIMEphoton<900) 

90 DIM $igmanu<900) 

100 DIM Ne < 900 ) 

110 REM 

120 REM 

130 REM F •=. c = Fine Structure Constant 

140 Fsc=l/137. 03604 

150 Hbarc=197. 32858 

160 Mnc sq = 938 . 95 

170 Mneutror. = 939. 5731 

180 Mprot on =9 38. 2796 

190 Arnu = 931 . 5016 

200 Mstar=.7*Mncsq 

210 J = 36 . 8 

220 Q= 1 7 

230 Eps i 1 on= . 0768 

240 INPUT "ENTER Z OF TARGET " , Zt 

250 INPUT "ENTER A OF TARGET", At 

260 Nt=At-Zt 

270 INPUT "ENTER Z OF PRO JECT I LE " , Zp 

280 INPUT "ENTER A OF PRO JECT I LE " , Ap 

290 Np=Ap-Zp 

300 IF Ap< 50 THEN Width=10.0 

310 IF Ap > = 50 THEN Width=4.5 

320 Frac 1=1. 95*EXP < - . 075*Zp ) 

330 F r ac 2 = Z p / A p 

340 IF Frac 1 < = Fr ac 2 THEN Fracprot on = Frac 1 

350 IF Frac2<Fracl THEN Frac prot on=Frac 2 

360 R 1 0t = 1 . 18*At^<l/3)+. 75 

370 R 1 0p= 1 . 1 8*Ap A < 1 /3 > + . 75 

380 Dee=0 

390 Bmin=R10t+R10p-Dee 

400 INPUT "WHAT IS KE/N OF PROJECTILE (MeV/N) ?",Tlab 

410 Gamma= 1 +T 1 ab/Mnc sq 

420 Vel =SQR < 1 - 1 /Gamma^-Z ) 

430 REM Gamma IS THE RELATIVISTIC GAMMA FACTOR OF PROJ 

440 REM Vel IS VELOCITY OF PROJ IN UNITS OF C (RELATIVISTIC BETA FACTOR) 

450 Si gmam=120*Np*Zp/(PI*Ap*Wi dth) 

460 Ro- 1 . 18*ApMl/3) 

470 U = 3*J*ApM-l/3)/Q 

480 Egdr = SGR ( 8 . 0* J*Hbarc - s 2 < Mst ar*Ro A 2 ) * 1 / < 1 +U- ( 1 +Eps i 1 on + 3 + U ) *Eps i 1 on/ ( 1 +Eps i 

1 on + U) ) ) 

490 REM 

500 REM NUMERICAL INTEGRATION OR PLOT 

510 REM 
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520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 * 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 


EphotorK 1 ) = . 001 
REM IF Et hreshgn>Et hreshgp 
Ephot onmax=70 
N p t s = 5 0 
REM 

REM Eint. is defined as the 
REM 


THEN Ephoton<l)=Ethreshgp 


integration or- plot in t e r v a 1 


Ei nt =< Ephot onmax-Ephot on < 1 > )x<Npts-l ! 
Sum = 0 
S u m p = 0 
Sumn=0 


REM 

REM 

REM 

FOR 


1=1 TO Npts 

Ephoton=Ephoton<l)+<I-l)*Eint 

Ephoton<I)=Ephoton 

S i gmanu=S i gmamx < 1 + < Ephot on /, -2-Egdr- A -2 ) A 2/' < Ephot on^2*W i dt > > 

Sigmanu(I)=Sigmanu 
Ecutoff=Hbarc*Gamma*Vel/Bmin 
G=EphotonxEcutoff 
CALL Bessel <G,K0,K1> 

N e = 2 * Z t 2 * F s c / < E p h o t o n * P I * V e 1 2 > * < G * K 0 * K 1 - . 5 * V e 1 x - 2 * G - A - 2 * C K 1 ■ s i 

Ne < I > =Ne 
Function=Si gmanu*Ne 
IF 1=1 THEN Func t i on= . 5*Func t i on 
IF I=Npts THEN Funct i on=. 5*Func t i on 
Sum=Sum+F unction 

Functionp=Fracproton*Function 
Funct ionn=< l-Fracproton)*Funct i on 
IF Ephoton<Ethreshgp THEN Func t i onp=0 
IF Ephot on<Ethreshgn THEN Functionn=0 
Sump=Sump+Func t i onp 
Sumn=Sumn+Func t i onn 


NEXT I 
REM 
REM 
REM 


IF 


Integralp=Eint*Sump 
Integraln=Eint*Sumn 
Integral=Integralp+Iritegraln 
Ephot onmax-Egdr< 40 THEN PRINT "WARNING: 
PRINT 


i ncre as e Ephotonm ax 1 


950 

PRINT 




960 

PRINT 

"Width <MeV>", Width 



970 

PRINT 

"Zt ",Zt 



980 

PRINT 

"At", At 



990 

PRINT 

"Zp",Zp 



1000 

PRINT 

"Ap" , Ap 



1010 

PRINT 

"KE/N < MeVxN ) " , T 1 ab 



1020 

PRINT 




1030 

PRINT 




1040 

PRINT 




1050 

PRINT 

"Lower limit of integration 

< MeV > 

" , Ephot on < 1 ) 

1060 

PRINT 

"Upper limit of integration 

< MeV > 

" , E p h o t o n m ax 

1070 

PRINT 

"Number of integration intervals 

is", Npts 

1080 

PRINT 

"Value of integration interval width <MeV>",Eint 

1090 

PRINT 




1100 

PRINT 




1 110 

PRINT 

" S i g m an u ( m b > " , S i g m an u 



1120 

PRINT 

" S i g m am < m b > " , S i g rn am 



1130 

PRINT 

"Ro <fm)",Ro 




2 - K 0 ■■■■■ 2 > ) 


ORIGINAL PAGE 15 

OP POOR QUALITY 
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1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

STOP 

END 


"IT , U 

" GDR Energy < MeV > n , Egdr 


"PRO J VELOCITY OBeta f ac t or > -un i t s of c",Vel 
"RELRTIVISTIC GAMMA FACTOR OF PROJ (MeV^N) Gamma 
"Ecutoff <MeV>", Ecutoff 

“10 percent charge radius of target <frn) ",R10t 
"10 percent charge radius of projectile <fm)",R10p 
"Dee", Dee 
" Bmi n <f m> " , Bmi n 


'COULOMB DISSOCIATION CROSS SECTION CSigmaww) (mb)", 


" S i gnia( gamma, p > 
"Si gma( gamma, n> 


< m b > " , I n t e g r a 1 p 
< mb > ", I nt egral n 


SUB Bessel <G,K0,K1> 
Al=3. 5156229 
A2=3. 0899424 
A3=l. 2067492 
A4=. 2659732 
A5=. 0360768 
A6=. 0045813 
A7=. 39894228 
A8=. 01328592 
A9=. 00225319 
A 1 0= . 00157565 
All=. 00916281 
A12=. 02057706 
R13=. 02635537 
A14=. 01647633 
A15=. 00392377 
A16=. 87890594 
A17=. 51498869 
A 1 8= . 15084934 
A19=. 02658733 
A20=. 00301532 
A21=. 00032411 
A22=. 39894228 
A23=. 03988024 
A24=. 00362018 
A25= . 00163801 
A26=. 01031555 
A27=. 02282967 
A28=. 02895312 
A29=. 01787654 
A30=. 00420059 
Bl=. 57721566 
B2=. 42278420 
B3=. 2306975 6 
B4=. 0348859 
B5=. 00262698 
B6=. 0001075 
B7=. 0000074 
B8= 1 . 25331414 
B9=. 07832358 
B10=. 02189568 


ORIGINAL) PAGE IS 
OF POOR QUALITY 


I nt egral 
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1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

I860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 


B1 1 = . 01062446 
B12=. 00587872 
B13=. 00251540 
B14=. 00053208 
B 1 5= . 15443144 
B16=. 67278579 
B17«. 18156897 
B18=. 01919402 
B19=. 001 10404 
B20=. 00004686 
B 2 1 = 1 . 25331414 
B22=. 23498619 
B23=. 03655620 
B24=. 01504268 
B25=. 00780353 
B26=. 00325614 
B27=. 00068245 
T=G/3.75 


ORIGINAL PAGE IS 
OF POOR QUALITY] 


IF G < - 3 . 75 THEN 1 0= 1 +fl 1 *T A 2+fi2*TM+fl3*T''-6+R4*T^8 + fi5*T^ 1 0 + fi6*T"' 1 2 

IF G >3 . 75 THEN I 0= 1 /SQR ( G > *EXP C G > * < fl7 + fl8^T + fl9^T A 2-fi 1 0.'T~3 + R 1 1 ^TM-fi 1 2^'T'-5 + 


fll3/T^6-fll4/T^7+fll5/T-8> 

I960 IF G< =3 . 75 THEN I 1 =G* < . 5 + fl 1 6*T-2+fl 1 7*T-4 + fl 1 S*T-"-6 + fl 1 9*T A 8+fl20*T A 1 0 + R2 1 *T~ 1 2 

> 

1970 IF G >3 . 75 THEN Il = l/'SQR<G)*EXP<G)*<fl22-fl23-'T-fl24/T A 2+fl25^T A 3-fl26^T' s 4 + fl27^T 
' s 5-fi28^T-*-6 + fl29/T' s 7-fl30^T' s 8> 


1980 S=G/2 

1990 IF G< =2 THEN K0 = -LOG < S > * I 0-B 1 + B2*S~2 + B3*SM + B4*S A 6 + B5*S A 8 + B6*S A 1 0+B7*S A 1 2 
2000 IF G >2 THEN K0= 1 /SQR < G > *EXP < -G > * < B8-B9/S+B 1 0/S^2-B 1 1 /S-3 + B 1 2/SM-B 1 3/S A 5 + B 
14/S^6) 

2010 IF G< =2 THEN K 1 =LOG < S > * 1 1 + 1 /G* ( 1 +B 1 5*S A 2-B 1 6*S~4-B 1 7*S-6-B 1 8*S~8-B 1 9*S A 1 0- 
B20*S A 1 2 ') 

2020 IF G >2 THEN K1 = 1/'SQR<G>*EXP<-G>*<B21 + B22/S-B23/S''2+B24/S A 3-B25/SM + B26^S A 5 
-B27/S-'6 ) 

2030 SUBEND 
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Sample Output 


Width ( MeV ) 

Zt 

fit 

Zp 

Rp 

KE/N (MeV/N) 


4.50 

82.00 

208.00 

26.00 

56.00 

1880.00 


Lower limit of integration ( M e V > .00 

Upper limit of integration (MeV) 70.00 

Number of integration intervals is 50.00 

Value of integration interval width (MeV) 1.43 


Sigmanu (mb) .56 

Sigmam (mb) 118.23 

Ro <fm) 4.51 

U 1.70 

GBR Energy (MeV) 18.40 


PRO-J VELOCITY (=Beta factor)-units of c .94 

RELRT I V I ST I C GRMMfl FRCTOR OF PROJ (MeV/N) 3.00 

Ecutoff (MeV) 42.95 

10 percent charge radius of target (fm) 7.74 

10 percent charge radius of projectile (fm) 5.26 

Dee 0.00 

Bmin (fm) 13.01 


COULOMB DISSOCIATION CROSS SECTION (Sigmaww) (mb) 1020.04 

S i gma( gamma, p ) (mb) 283.00 

S i gma( gamma, n ) (mb) 737.05 
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Table 1. Resonance Widths and Proton Branching Ratios 


Nucleus 

r, MeV 

9p 

Reference 3 

Present work 

Reference 3 

Present work 

7 Li 


10.0 


0.43 

9 Be 


10.0 


0.44 

12 C 

8.0 

10.0 

0.5 

0.5 

16 0 

10.0 

10.0 

0.5 

0.5 

18q 

12.0 

10.0 

0.4 

0.44 

20 Ne 

10.0 

10.0 

0.5 

0.5 

28 Si 

10.0 

10.0 

0.5 

0.5 

32 S 


10.0 


0.5 

40 Ar 

10.0 

10.0 

0.45 

0.45 

40 Ca 

10.0 

10.0 

0.5 

0.44 

48 Ti 


10.0 


0.37 

54 Fe 

3.0 

4.5 

0.7 

0.28 

56 Fe 

5.0 

4.5 

0.28 

0.28 

58 Ni 

10.0 

4.5 

0.5 

0.24 

63 Cu 

5.0 

4.5 

0.28 

0.22 

9° Zr 

4.0 

4.5 

0.05 

0.10 

107 Ag 

5.0 

4.5 

0 

0.06 

106 Gd 

4.0 

4.5 

0 

0.02 

181 Ta 


4.5 


0.01 

197 Au 

3.5 

4.5 

0 

0.01 

208 Pb 

3.9 

4.5 

0 

0 

238JJ 

5.0 

4.5 

0 

0 



Table 2. Giant Dipole Resonance Particle Thresholds 



Proton threshold, 
MeV 

Neutron threshold, 
MeV 

Nucleus 

Reference 3 

Present work 

Reference 3 

Present work 

12 c 

15.46 

0.001 

18.74 


160 

11.62 

0.001 

15.67 


18 0 

15.44 

0.001 

8.05 


40 Ar 

12.02 

0.001 

9.87 


56 Fe 

9.67 

0.001 

11.20 


197 Au 

5.27 

0.001 

8.07 






Table 4. Calculated Total Electromagnetic Absorption Cross Section for 
1.88 GeV/N 56 Fe Incident Upon Various Targets 


Projectile 

Energy, 

GeV/N 

Target 

^EM(W), mb 

(a) 

OEM* m b, for d = 0 fm 

Reference 3 

Present work 

56 Fe 

1.88 

l Li 

2 

1.9 

2.1 



9 4 Be 

3 

3.3 

3.7 



‘ic 

7 

7.3 

8.0 



32 q 
16 s 

46 

46 

52 




130 

140 

156 



X 4?Ag 

306 


377 



^Ta 

629 

717 

830 



2 08 Pb 

834 

901 

1020 



2 i|u 

1008 

1105 

1250 


°This column represents the isotope-averaged calculations of Westfall et al. (ref. 6). 



Table 5. Calculated Total Electromagnetic Reaction Cross Sections for 12 C and 
16 0 Incident Upon Various Targets 


Projectile 

Energy, 

GeV/N 

Target 

* Final state 

<t E m(HL), mb 
(a) 

ctem ? nib, for d = 0 fm 

Reference 3 

Present work 

12 C 

2.1 

2°8pb 

u C + n 

50 ± 18 

54 

68 




u B + p 

50 ±25 

60 

68 


1.05 


n C + n 

38 ±24 

32 

43 




n B+p 

50 ±26 

36 

43 

16 o 

2.1 


15 0 + n 

50 ±25 

78 

99 




15 N + p 

97 ± 17 

87 

26 

12 c 

2.1 

108 Ag 

n C + n 

22 ± 12 

21 

26 




n B + p 

20 ±12 

23 

26 


1.05 


u C + n 

22 ± 12 

13 

17 




u B + p 

25 ±20 

15 

17 

16 q 

2.1 


15 0 + n 

26 ± 13 

30 

37 




15 N + p 

29 ± 18 

33 

37 

12 C 

2.1 

64 Cu 

u C + n 

10 ±6 

9 

11 




U B + p 

4 ± 8 

10 

11 


1.05 


n C + n 

10 ±7 

5.9 

7.4 



1 

n B + p 

5 ± 8 1 

6.5 

7.4 

16 o 

2.1 


15 0 + n 

10 ±7 

12.7 

16 




15 N + p 

14 ±9 

14 

16 

12 C 

2.1 

to 

> 

“C + n 

0 ± 3 

2.1 

2.5 




n B + p 

0 ± 3 

2.3 

2.5 


1.05 


n C + n 

1 ± 3 

1.5 

1.8 




11 B + p 

1 ± 3 

1.6 

1.8 

16 o 

2.1 


15 0 + n 

0 ± 3 

2.9 

3.6 




15 N+p 

0 ± 0 

3.2 

3.6 

12 C 

2.1 

12 c 

n C + n 

0 ± 1 

0.50 

0.58 




11 B + p 

0 ± 3 

0.54 

0.58 


1.05 


n C + n 

0 ± 2 

0.36 

0.43 




n B + p 

0±1 

0.40 

0.43 

16 o 

2.1 


15 0 + n 

0 ± 2 

0.70 

0.83 




15 N + p 

0 ± 3 

0.76 

0.83 


a This column represents the measurements (isotope averaged) of Heckman and Lindstrom (ref. 7). 
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Table 6. Calculated Total Electromagnetic Reaction Cross Sections for 18 O at 
1.7 GeV/N Incident Upon Various Targets 







(TEm, mb, for d = 0 fm 

Projectile 

Energy, 

GeV/N 

Target 

Final state 

mb 

(«) 

Reference 3 

(9p = 0-2) 

Present work 

18 0 

1.7 

48 Ti 

17 0 + n 

8.7 ±2.7 

12 

ii 




17 N + p 

0.5 ±1.0 

2 

9 



208 Pb 

17 0 + n 

136 ±2.9 

123 

112 




17 N + p 

20.2 ±1.8 

24 

89 


; 

238^ 

17 0 + n 

140.8 ±4.1 

151 

136 



I 

17 N + p 

2.51 ± 1.6 

30 

109 


“This column represents the measurements (isotope averaged) of Olson et al. (ref. 8). 


Table 7. Target Fragmentation — Calculated Total Electromagnetic Reaction Cross Sections 

for Various Projectiles Incident Upon 197 Au 







cteM) mb, for d = 0 fm 

Projectile 

Energy, 

GeV/N 

Target 

Final state 

<t em (M), mb 
(a) 

Reference 3 

Present work 

12 C 

2.1 

197 Au 

196 Au + n 

66 ±20 

37 

39 

20 Ne 

2.1 



136 ± 21 

97 

104 

40 Ar 

1.8 



420 ± 120 

278 

299 

56 Fe 

1.7 



680 ± 160 

546 

588 


“This column represents the data of Mercier et al. (ref. 9). 
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reaction 
cross section, 
mb 



Figure 1. Theoretical, parameterized, and experimental photoneutron reaction cross sections for 12 C. 
Theoretical T = 8 MeV; Parameterized T — 10 MeV; Theoretical g v = Parameterized g p = 0.5; thus, 
photoneutron and photoproton cross sections are identical. 



Figure 2. Theoretical, parameterized, and experimental photoneutron reaction cross sections for 16 0. 
T = 10 MeV; Theoretical g v — Parameterized g v = 0.5; thus, the photoneutron and photoproton cross 
sections are identical. 
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Virtual photon energy, MeV 

Figure 3. Theoretical, parameterized, and experimental photoneutron reaction cross sections for 28 Si. 
T = 10 MeV; Theoretical g p — Parameterized g p = 0.5; thus, the photoneutron and photoproton cross 
sections are identical. 



Virtual photon energy, MeV 

Figure 4. Theoretical, parameterized, and experimental photoneutron reaction cross sections for 58 Ni. Theoret- 
ical T = 10 MeV; Parameterized F = 4.5 MeV; Theoretical g p = 0.5; thus, the theoretical photoneutron and 
photoproton cross sections are identical; Parameterized g p = 0.24. Resultant parameterized photoproton 
cross section is also displayed. 
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Figure 5. Theoretical, parameterized, and experimental photoneutron reaction cross sections for ^Zr. 
Theoretical T = 4 MeV; Parameterized T = 4.5 MeV; Theoretical g p = 0.05; Parameterized g p = 0.1. 



Figure 6. Theoretical, parameterized, and experimental photoneutron reaction cross sections for 160 Gd. 
Theoretical T = 4 MeV; Parameterized T = 4.5 MeV; Theoretical g p = 0; Parameterized g p — 0.02. 
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Figure 7. Theoretical, parameterized, and experimental photoneutron reaction cross sections for 197 Au. 
Theoretical T = 3.5 MeV; Parameterized F = 4.5 MeV; Theoretical g p = 0; Parameterized g p — 0.01. 



Figure 8. Theoretical, parameterized, and experimental photoneutron reaction cross sections for 208 Pb. 
Theoretical T = 3.9 MeV; Parameterized F = 4.5 MeV; Theoretical g p = Parameterized g p = 0. 


20 



Figure 9. Theoretical, parameterized, and experimental photoneutron reaction cross sections for 238 U. 
Theoretical r = 5 MeV; Parameterized T — 4.5 MeV; Theoretical g p = Parameterized g p = 0. 


Photonuclear 
reaction 
cross section, 
mb 



Figure 10. Theoretical, parameterized, and experimental photoproton reaction cross sections for 18 0. 
Theoretical T — 12 MeV; Parameterized T = 10 MeV; Theoretical g p — 0.40; Parameterized g p = 0.44. 
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Photonuclear 
reaction 
cross section, 
mb 



Virtual photon energy, MeV 

Figure 11. Theoretical, parameterized, and experimental photoneutron reaction cross sections for 18 0. 
Theoretical T = 12 MeV; Parameterized T = 10 MeV; Theoretical g n = 0.60; Parameterized g n = 0.56. 



Figure 12. Theoretical, parameterized, and experimental photoproton reaction cross sections for 54 Fe. 
Theoretical T = 3 MeV; Parameterized Y = 4.5 MeV; Theoretical g p = 0.70; Parameterized g p = 0.28. 


22 




Figure 13. Theoretical, parameterized, and experimental photoneutron reaction cross sections for 54 Fe. 
Theoretical T = 3 MeV; Parameterized V = 4.5 MeV; Theoretical g n — 0.30; Parameterized g n = 0.72. 
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